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by 
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SUMMARY 


In addition to the forces that influence the motion of the center of mass of a satellite, 
various forces are also present which influence the motion about the center of mass. 
These torque -producing forces may result in a change in satellite orientation that affects 
the thermal balance, operation of solar cells, various scientific measurements, and 
drag forces. This latter effect results in a coupling between the equations describing 
the orientation with those describing the orbital motion. Therefore, strictly speaking, 
the orientation problem is part of the orbital problem. 

To solve the orientation equations it is necessary to determine the nature of the 
torque-producing forces* This is done by observing the change in orientation of several 
Explorer satellites where the drag force is sufficiently small so the orbital dependence 
on orientation is negligible. It was found that permanent magnetic moments in the sat- 
ellites were the dominant effect* responsible for the observed changes in orientation. 
Gravitational torques are also significant. The changes in orientation of the satellites 
considered were explained extremely well by these two effects, hence it is concluded 
that other effects are not significant for similar satellites. 

Various approximations used in this study greatly reduce the effort required to 
integrate the orientation equations and do not require the simultaneous solution of the 
orbital equations. The findings, however, are applicable to the simultaneous solution 
of both the orbital and orientation sets of equations. 


SECTION I. INTRODUCTION 


The most general treatment of satellite motion requires six second-order differen- 
tial equations to describe the motion. Three equations describe the motion of the center 



of mass in phase space, and the remaining three equations describe the orientation and 
rates of three body-fixed axes relative to a space-fixed system. In orbital calculations 
it is usually assumed that this set of six equations separates into two sets with the orbital 
set independent of the set describing the orientation. This is not rigorously true if the 
satellite in question has any asymmetry, since the drag forces will depend on the orien- 
tation of the satellite. 

It is desirable for several reasons to formulate the full six degree of freedom prob- 
lem and integrate the equations describing the orientation along with the orbital equations. 
It may be argued that since the previously mentioned drag dependence is usually small 
compared to the uncertainties in atmospheric density, such an exercise is largely aca- 
demic. While this is a valid argument there are many satellite experiments in which it 
is just as important, or even more important, for the experimenter to have knowledge 
of the satellite orientation as it is to know the position of the center of mass. This is 
particularly true in the case of the Explorer XI gamma ray astronomy satellite. 

In most cases where it is necessary to have orientation information, the approach 
has been to derive such information from a combination of solar and earth horizon sen- 
sors contained in the satellite, or by analyzing the recorded radio signal strength pat- 
terns. The latter has been particularly successful with Explorers IV, VII, and XI. 

While it has been possible to determine satellite orientation a posteriori, it would 
be desirable to be able to predict the satellite orientation a priori. Since the external 
torque-producing forces are coupled closely to the motion of the center of mass, and 
since the motion of the center of mass depends to a lesser degree on the orientation, it 
would seem reasonable to consider the entire six degree of freedom formulation in the 
integration of the equations of motion. A convenient formulation of the full six degree 
problem has been suggested by Lundquist and Naumann [Ref. 1] in which the orientation 
of the body-fixed axes are expressed in terms of the Euler quaternion parameters 
[Ref. 2]. This formulation has several advantages over the more conventional Euler 
angle formulation in that it yields symmetrical equations in which no trigonometric 
functions are involved. These parameters also have certain advantages over the Caley- 
Klein parameters [Ref. 3] in that complex numbers are not involved. The disadvantage 
is that they are not independent and exceed by one the number of degrees of freedom. 

This difficulty is alleviated by the use of a Lagrange undetermined multiplier and the 
constraint relationship that the sum of the squares of the parameters is unity. Several 
computer programs based on this formulation have been developed by the Computation 
Division at the Marshall Space Flight Center and by Cunningham at the University of 
California. 

In order to use such programs to predict satellite motion about the center of mass, 
it is necessary to know the nature of the external torques acting on the satellite. Anal- 
ysis of the motion of several of the Explorer satellites has yielded some knowledge con- 
cerning the nature of these torques. The purpose of this paper is to discuss some of the 
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experimental satellite orientation data, and attempt to explain the nature and origin of 
the torques responsible for the observed motion about the center of mass. 

The satellites to be considered are Explorer XI, Explorer IV, Explorer VIII and 
Explorer VII. This particular order was chosen because Explorer XI has yielded the 
most accurate and extensive orientation information.. Sufficient information has been 
obtained to allow the extraction of the torques by numerical differentiation, hence, the 
determination of significant effects. Much less orientation data are available for 
Explorer IV, but in the light of the Explorer XI analysis, its motion is now fairly well 
understood. Very little data exist on Explorers VIII and VII, but again using results 
from the previous analyses a motion is derived that is consistent with known data. 


SECTION II. COMPUTATIONAL TECHNIQUES 


Before discussing the various satellites individually, it would be well to discuss 
various sources of torques and the various approximations used in computing their 
effects. 


To begin, it will be assumed that all of the satellites have cylindrical symmetry 
and are initially set spinning about their axes of symmetry. The moment of inertia 
about the axis of symmetry is denoted as I 3 . . If a transverse axis has a moment of 
inertia I 2 less than I 3 , the satellite is said to be stable and will remain spinning about 
I 3 . If I 2 is greater than I 3 , as is the case in Explorer XI and Explorer IV, the satellite 
is in a maximum rotational energy configuration and is unstable. Due to internal 
energy dissipation and angular momentum conservation, the satellite longitudinal axis 
will precess about the angular momentum vector in a cone. The half-angle of this 
precession cone a is related to the longitudinal spin <n 3 by 


cos a 


I 3 

L 


(1) 


where L is the total angular momentum. As energy is dissipated, w 3 -*■ 0, a — n/2, 
and the satellite approaches the minimum energy configuration in which the rotation 
is a tumbling motion with the axis of symmetry rotating in a plane perpendicular to the 
angular momentum vector. 

A. GRAVITATIONAL GRADIENT TORQUES 

The gravitational torque on a satellite may be found by expanding the potential 
energy about the center of mass in a Taylor series and differentiating with respect to 
the generalized angles expressing the orientation of the body-fixed axes. The result 
given by Roberson [Ref. 4] and others [Ref. 5] is 
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( 2 ) 


mk 2 a a a a 

L g = 3-p- (Ii - I 3 ) ( A • R) (AxR) 

where mk 2 is the gravitational constant, r is the radius vector magnitude, A is a unit 
vector along the axis of cylindrical symmetry, and R is a unit vector along the radius 
vector. 

Assuming a half-cone angle a, this torque may be time averaged over a pre- 
cessional or tumble cycle while holding R constant. This results in 



3_ mk 2 
2 'r 3 ” 


(Ii - I 3 ) (R ’ L) (RxL) (1-3 cos 2 a) 


(3) 


A 

where L is a unit vector directed along the angular momentum vector. 

Assuming a circular orbit, this L„ is time averaged over an orbital revolution 

A (j A 

while holding L constant. This is done by expressing R in terms of orbital inclination 
i, right ascension of the ascending node 0 , and mean anomaly M. 


/ cos M cos 0 - cos i sin 0 sin M 

A 

R = cos M sin <p + cos i cos 0 sin M 
y sin i sin M 

The unit vector along the orbital angular momentum is 

( sin i sin 0 \ 

- sin i cos 0 
cos i 


A 

12 = 


\ 


/ 


(4) 


(5) 


The time average is found by 

.2 


2 7T 


l g - 


47T 


mk a c /* A A A A 

^ (I 2 - I3) (1-3 cos 2 a) J (R: L) (Ex L) dM . 
r 0 


( 6 ) 


A 

Expressing R in terms of 0 , i, and M from equation 4 , integrating, and ex- 
pressing the resulting combinations of 0 and i in terms of components of 12 using 
equation 5, results in 
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(7) 


= k q (L-Q) 


A. 

(Lx!l) 


( i - 3 cos 2 o;) 


where 


k g = T ~h)- (8) 

B. MAGNETIC TORQUES 

If a satellite has a permanent magnetic moment projection, M^, on the axis 
of rotation, an interaction torque will arise from the geomagnetic field. This torque 
is given by 


L„ = M t x B = M t (L x B) . 
ML L 


(9) 


A component of magnetic moment along the axis of symmetry, M 3 , results in 
M = M 3 cos a. Due to rotational symmetry, a transverse component of magnetic 
moment does not contribute unless co 3 = JD. If this happens, the M^ is M t or the pro- 
jection of the transverse component on L. 

To facilitate computation, an earth-centered dipole representation of the geo- 
magnetic field is employed, 


B 


/ - 3R x R z\ 



-3 

1 


R R 

y z 

- 3R 2 
z 
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( 10 ) 


where B^ is the scalar field at the equator at 
This is given by 


r km from the center of the earth. 


B 


N 



(ii) 


where y is the geomagnetic moment taken as 8. 1 x 10 55 cgs units. 

Assuming a circular orbit, the magnetic torque is time averaged over an orbital 
revolution in the same manner as before, giving 


V 
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The constant M T B__ 
L N 


is termed the magnetic couple. 


(12) 


C. INDUCED MAGNETIC TORQUES 

For a satellite having a long cylindrical geometry, it is possible that a mag- 
netic moment may be induced that alternates in sign during a tumble cycle [Ref. 6]. 
Torques produced by this process would tend to add rather than cancel because of 
rotational symmetry. The induced moment is 

M = -^rr ~ - U (B -A) A ( rationalized MKS units) (13) 

I ro 

where u is the relative permeability; u is the permeability of free space; and U 
r o 

is the volume of material in the walls of the cylinder. Let 

K = — - U, (14) 

1 "o 


Kj (A •' B) (A x B) . (15) 


then the torque is 


Time averaging over a tumble cycle results in 

Lj = - y Kj (L-B) (L XB) . (16) 

Time averaging over an orbital revolution, as before, results in the expressions 
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L I = 




(17) 


where 


where 


G i * B ij L j ; 1 = 15 2 > 3; j = 1 > 2 ’ 3 


Li — L * Lo — L j Lo — L 
1 x y d z 


and 



X y 


ft ft (1 - 3 ft 2 ) 
x y z 


fi S2 (1-3S2 2 ) n £2 (4-9/3 2 ) 

8 x y v z ; 8 x z 


9 (« 2 + 3£2 2 fi 2 ) ■— (4 - 9 /3 2 ) 


8 ' x 
3 


y z 


8 y z 


£2 x ^ z (4- 9^ 2 ) -J- £2 y 5^ (4- 9/3 2 ) 1 - 3 )3 2 + y j£ 


/3 2 = n 2 + n 2 

x y 


D. DRAG TORQUES 

The drag force on a satellite is given by 


D 


2 °D PSV V ’ 


(18) 


where is the drag coefficient usually taken to be 2. 0, p is the atmospheric density, 
and S is the area projected normal to the velocity vector V . The drag torque is 


l d = a x f d 


(19) 


where a is the displacement of the center of pressure from the center of mass. 
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For symmetrical satellites, such as Explorers VII and VIH, the a is for all 
practical purposes zero and the drag torque is not significant. For'asymmetrical sat- 
ellites, such as Explorer XI, the drag torque is significant only before the satellite 
has reached the tumbling mode. The first order effects cancel by rotational symmetry 
after a becomes tt/2. 

For an eccentric orbit, such as Explorer XI, the atmospheric density is signifi- 
cant only near perigee. Hence, the V is taken as the perigee velocity vector. The 
time average torque over an orbital revolution is found by 

L d = - K d (Lx V p ) (20) 

where 

k d =i c D sv‘.ap 


and p is the time average atmospheric density given by 

2 7T 

p = p ( r) d M • 

0 

E. ACCURACY OF APPROXIMATIONS 

In order to determine the errors introduced in assuming circular orbits and 
the earth-centered dipole representation for the geomagnetic field, numerical integra- 
tions were carried out over actual orbits of Explorer XI, using the Finch and Leaton 
48-term Legendre Polynomial representation for the magnetic field [Ref. 7], 

It was found that the gravitational torques computed using the approximate methods 
differed from those computed by integration over an entire orbit by a maximum of 1%. 

The permanent magnetic torques differed by approximately 10% using a value of 
B of . 191 gauss computed from equation 11. However, it was noted that the values 
obtained from the numerical integration were always larger than those from the 
approximate method. Assuming a of . 205 gauss reduced this difference to 
approximately ±3%. This is not serious in this analysis since the is not known 
accurately and the product must be assumed. 

The computation of the induced torques suffers somewhat more from the approxi- 
mations. The maximum observed deviation was 14% of the maximum value, but 


8 



generally the deviations were no more than ±8% of the maximum value. The spread 
was such that no apparent improvement could be gained by slightly altering the value 

° f V 


SECTION III. EXPLORER XI 


The orientation history of Explorer XI was obtained by determining the angle a 
line from the tracking station to the satellite makes with angular momentum vector 
by analyzing the recorded radio signal strength variations [Ref. 8]. Many such angles 
with their respective line of sight vectors can be determined during a set of successive 
passes. These data result in an overdetermined set of non-linear equations that are 
solved by an iterative process to yield the angular momentum unit vector that best 
fits the observed data in the sense of least squares. The observed orientation is 
plotted in terms of right ascension and declination in Figure 1. 

Beginning at day 35 after 27. 0 April 1961, these plots are fitted in segments by 
polynomials which are then differentiated with respect to time to yield the torques 
responsible for the motion. 

Since the gravitational torques can be computed accurately from the known moments 
of inertia of the satellite (Table 1) , they are subtracted from the observed torques and 
the results are plotted in Figure 2. This resulting torque must then be caused by 
either a permanent or induced magnetic moment. The induced torques are computed 
from equation 12 assuming a value of -384 dyne-cm for -g- B-^- 2 which was chosen 

to yield the approximate observed amplitude. It may be seen from Figure 2 that 
the induced torques have the same general shape as the required torques, but the 
correlation is not very convincing. A further difficulty arises in the negative sign that 
must be assumed for since should be positive. Several points are shown that 
were computed from the detailed integration procedure. It is obvious that the errors 
introduced in the approximations are small compared to the differences between the 
induced torques and the required torques. 

On the other hand, the points resulting from an assumed value of 148. 14 dyne-cm 
for the permanent magnetic couple fall right on the plots of the required torques. 

Therefore, it is concluded that a permanent transverse magnetic moment is the dom- 
inant effect responsible for the observed motion. 

The observed torques together with the gravitational and permanent magnetic 
torques are plotted in Figure 3 for day 35 through day 180. It is seen that the sum 
of gravitational and magnetic torques equals almost exactly the observed torques, the 
residuals being generally less than five dyne-cm. It is believed that the errors in 


9 



RIGHT ASCENSION (degrees) DECLINATION (degrees) 



DAYS AFTER APRIl ?70 UT IQfil 

FIGURE i. OBSERVED ORIENTATION OF THE ANGULAR MOMENTUM VECTOR OF EXPLORER XI 


TABLE 1 


DYNAMICAL PROPERTIES 


QUANTITY 

EXPLORER XI 

EXPLORER IV 

EXPLORER VII 

EXPLORER VIH 

Time Origin 

27. 0 Apr 1961 

26. 0 Jul 1958 

13. 0 Oct 1959 

3. 0 Nov 1960 

Ii (gm-cm 2 ) 

1.627 x 10 s 

4.752 x 10 7 

1.08 x 10 7 

2. 162 x 10 7 

I 3 (gm-cm 2 ) 

0. 040 x 10 8 

0. 056 x 10 7 

2.93 x 10 7 

2.936 x 10 7 

Ii-I 3 (gm-cm 2 ) 

1.587 x 10 8 

4. 696 x 10 7 

-1.85 x 10 7 

-0. 774 x 10 7 

r (km) 

7512 

7616 

7191 

6964 

K^(dyne-cm) 

111.92 

31.78 

-14. 87 

-6.85 

M 1 B^(dyne-cm) 

148. 14 

N/A 

N/A 

N/A 

M 3 B N (dyne-cm) 

450 

-3500 

75 

-200 

Mj (amp-m 2 ) 

.7756 

N/A 

N/A 

N/A 

M 3 (amp-m 2 ) 

2. 35 

-19. 125 

.345 

-.837 

0 O (deg) 

253.912 

65. 81 

149.217 

16. 356 

$> (deg/ day) 

-5. 0036 

-3.6505 

-4. 1745 

-3. 3698 

i (deg) 

28.8 

51.0 

50, .27 

49.98 

L (gm cm 2 /sec) 

1.276 x 10 8 

4. 568 x 1'0 7 

1. 104 x 10® 

3. 073 x 10 8 

AL/ day 

-0. 223 x 10 7 * ** 

-0. 0076 x 10 7 

N/A 

-0. 032 x 10 7 


-0. 0143 x 10 7 


* for first 20 days before tumbling commenced 

** after 26 days when a became ir/2 



L% (DYNE-CM) L\ (DYNE-CM) L # x DYNE-CM 



FIGURE 2. ANALYSIS OF THE PERMANENT AND INDUCED MAGNETIC 


TORQUES COMPARED TO THE DIFFERENCE BETWEEN THE 
OBSERVED TORQUE AND GRAVITATIONAL TORQUE FOR 
EXPLORER XI 
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obtaining the torques from the differentiation of the fitted observations are primarily 
responsible for these deviations, since the greatest deviations usually occur near the 
ends of the fitted segments. 

It is known from solar cell output data that Explorer XI had no o> 3 after day 26. 
This probably occurred because one transverse axis has a slightly greater moment of 
inertia; hence, there is a preferred axis of tumble! Also, a mercury toroidal damper 
near the tail is quite effective in removing any residual roll about the longitudinal axis. 
Therefore, it ijs quite likely that a transverse magnetic moment can have a constant 
projection on L. 

There is al&o evidence that a transverse magnetic moment did indeed exist in the 
attached stainless steel rocket case. Handling difficulties with the loaded rocket 
casing prevented an accurate measurement of the magnetic moments, but a small hand 
compass held near the case showed evidence of a transverse as well as a longitudinal 
magnetic moment. 

Although it has been shown that the permanent magnetic moment is the dominant 
effect, it is still necessary to estimate the magnitude of the induced effect. This is 
extremely difficult to do in a meaningful manner, since the relative permeability is not 
well known. In ferrous materials, the permeability depends strongly on the applied 
field as Well as on the past magnetic history of the material. Taking the volume of the 
ferrous material in the walls of the rocket casing to be 0. 470xl0 -3 m 3 (the payload 
portion is fabricated from aluminum) with a permeability of 100 in a field of 0. 2 x 10 -4 
weber/m 2 results in a K B 2 of 150 dyne-cm. However, the value of 100 for ju , 
which is generally considered the initial permeability of 410 stainless steel applies for 
a flux density of 200 gauss in the material. It is probably considerably less for an 
applied field of ±0. 2 oersted experienced in the geomagnetic field. There may well be 
a small component of induced torque, but this would act in the same general manner as 
the torque due to the permanent magnetism, but with opposite sign. This would prob- 
ably be compensated for by a slightly greater M than was assumed. 

Since it appears that the gravitational and permanent magnetic torques are mainly 
responsible for the motion, an attempt was made to supplement the few orientation data 
available between launch and day 20 before the satellite began tumbling. This was 
done by starting at the injection orientation and integrating the motion, assuming 
various values of M 3 to see if a motion that was consistent with the observed data could 
be made to result. It was found that a couple of M 3 B^. = 450 dyne-cm gave a reasonable 
fit of the observed orientations. The addition of a drag couple of 50 dyne-cm further 
improved the fit. The results are shown in Figure 4. 

Also, the motion was integrated from day 26 to day 90 and compared with observed 
data. As may be seen in Figure 5, there is excellent agreement between computed and 
observed results. 
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FIGURE 4. COMPUTED ORIENTATION OF THE ANGULAR MOMENTUM VECTOR COMPARED 
TO OBSERVATIONS FOR THE FIRST 20 DAYS OF EXPLORER XI 
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SECTION IV. EXPLORER IV 


The orientation of Explorer IV was deduced in much the same manner as Explorer 
XI, except that in Explorer IV this analysis was done as an afterthought. Therefore, 
the choice of tracking antenna and recording techniques was far from optimum for such 
an analysis. As a result, the amount and accuracy of the orientation data are much 
less than for Explorer XI. The data were sufficient to explain the temperature meas- 
urements in terms of area presented to the sun [Ref. 9], Also, the motion of a direc- 
tional scintillation counter relative to the magnetic field was deduced [Ref. 10] which 
allowed a determination of the directional flux of the trapped corpuscular radiation in 
the radiation belts [Ref. 11], The orientation measurements are shown in Figure 6. 

It appears that the orientation changes at approximately the same rate as observed in 
Explorer XI. However, the torques responsible for this motion remained unexplained 
until quite recently. The possibility of induced magnetic torques was suggested by 
Colombo [Ref. 12] , but starting at day 32. 0 after 26. 0 July 1958 and integrating the 
motion produced by gravitational and induced torques resulted in a motion that is not 
consistent with the observed data shown by the short dashed line in Figure 6. 

Integrating the torques produced by a permanent tranverse magnetic couple of -100 
dyne-cm, similar to that found in Explorer XI, produced a motion that apparently 
matches the observed data for several days before and after 32. 0. A serious difficulty 
arises with this assumption, however, since it was known that Explorer IV had a 
residual a> 3 which would tend to cancel any effect due to a transverse magnetic moment. 
Evidence for this roll is seen as modulations in the observed count rate of a directional 
scintillation detector whose axis is directed along a transverse axis [Ref. 13]. Fig- 
ure 7 shows the roll period history deduced from these count rate modulations. Note 
that the roll period increases rapidly after day 30. In Figure 6, the change in orien- 
tation seems to be much less after this time. From the roll period data, it can be 
determined that the half-cone angle oi increased rapidly the first few days, but ap- 
proached 7t/2 very slowly. In fact, from equation 1 it may be shown that oi is about 
83° on day 30. It then becomes clear that a longitudinal magnetic moment will give a 
constant projection on L that will produce the same effect as a transverse magnetic 
moment. A rather large value of -3500 dyne-cm must be assumed for this couple to 
give a projected component sufficient to account for the motion between day 25 and 
day 32, but this assumption appears to explain the observed data rather well. After 
day 32, the torques from gravitational effects appear to be sufficient to explain the 
motion in a reasonable manner. 

There may also be some significant second-order torques produced by the large 
longitudinal magnetic moment after day 32. These second-order torques arise from 
the fact that the geomagnetic field is changing because of geographical changes of the 
satellite during the 7-second tumble cycle. Such a change prevents the longitudinal 
magnetic torques from averaging to zero due to rotational symmetry. Using the 
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RIGHT ASCENSION (DEG) ' DECLINATION (DEG) 



FIGURE 6. ANALYSIS OF THE COMPUTED ORIENTATION OF 
EXPLORER IV DUE TO VARIOUS ASSUMPTIONS 
COMPARED TO OBSERVATION 




formulation of the full six degree of freedom problem discussed in the introduction, 
it was found that the net second-order effect was approximately 10“ 3 of the maximum 
instantaneous magnetic torque. This effect would therefore become significant only if 
there exists a large magnetic moment having rotational symmetry and if the other 
torque-producing effects are small. 

Some justification is necessary for the fact that a couple of -3500 dyne-cm corre- 
sponding to a longitudinal moment of 19 amp-m 2 is required to explain Explorer IV 
motion, whereas only 450 dyne-cm or 2. 35 amp-m 2 is required for Explorer XI. 

This can be resolved, qualitatively at least, by the fact that the instrumented portion 
of Explorer IV is stainless steel rather than aluminum. Furthermore, this portion 
was placed on a magnetically driven shake table for vibrational testing prior to flight. 
This quite probably resulted in a rather large permanent magnetization. Taking the 
volume of ferrous material to be 450 cm 3 , a flux density of approximately 500 gauss in 
the material is required to produce a magnetic moment of 19 amp-m 2 . The saturation 
of 410 stainless steel is on the order of 16 kilogauss and the retentivity is approx- 
imately 8 kilogauss for 10 kilogauss applied. Hence, 500 gauss is well within the 
magnetic capabilities of the material. 

It is curious that Explorer IV began tumbling very shortly after injection but re- 
quired 30 or more days to approach 7r/ 2. The longitudinal roll u> 3 apparently never 
reached zero and, in fact, even appears to increase on day 33. This increase is possi- 
bly due to an interaction torque from a transverse magnetic moment. On the other 
hand, Explorer XI had a very small half-cone angle until about day 20, at which time 
it opened rapidly and lost all its u> 3 within a few days. These facts may be explained 
by the consideration that Explorer IV did not have a mercury damper and the dissi- 
pative forces became less effective in removing energy as co 3 became small. It is 
difficult to understand why Explorer XI, with a mercury damper, took so long to begin 
tumbling. There is evidence that a tape recorder inside the satellite was accidentally 
commanded on and left running after launch. There is also evidence that this recorder 
failed on about the day the tumbling started. It has been suggested that the recorder, 
through bearing friction, actually put rotational energy into the satellite which over- 
came the energy dissipated by mechanical flexing [Ref. 14]. This presumably allowed 
the satellite to remain in its maximum energy configuration despite dissipative losses. 


SECTION V. EXPLORER VHI 


Unlike the satellites discussed previously, Explorer VHI was of a stable configura- 
tion. No attempt was made to analyze the orientation by radio signal strength tech- 
niques, but the satellite contained a crossed-slit solar aspect angle measuring device 
along with a horizon sensor. Presumably the data from both of these devices would 
determine the orientation of the spin axis. The solar aspect device yielded excellent 
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data, accurate to about 1 degree. However , the horizon scan data are apparently 
difficult to interpret. The only published orientation data are the angles between the 
spin axis and the sun for the first 30 days [Ref. 15], 

From the previous analysis, it was felt that the torques acting on Explorer VIII 
must be predominately due to a permanent magnetic moment along the spin axis and 
gravitational torques. Therefore, an attempt was made to integrate the motion from 
the injection condition to produce the observed solar aspect angle. 

The first measured solar aspect angle did not agree with an assumed initial spin 
axis coinciding with the injection velocity vector; hence, it was necessary to assume 
a spin axis injection condition a few degrees from the velocity vector at injection. This 
can be reconciled by a slight misalignment of the fourth stage velocity vector with the 
total injection velocity, or by assuming that the satellite was precessing slightly when 
the fourth stage was jettisoned. 

Various values of M 3 B^ were assumed until it was found that M 3 B^ = -200 dyne-cm 
.would produce a computed solar aspect angle shown in Figure 8. This computed result 
compares with the measured solar aspect angle well within the limits of accuracy. The 
resulting trace of the spin axis is shown in Figure 9. 


SECTION VI. EXPLORER VII 


A crude determination of Explorer VII orientation was accomplished by noting the 
times at which the fading caused by Faraday rotation of the 20 me/ sec signal strength 
became the most severe [Ref. 16]. The turnstile antenna emits omnidirectional 
radiation that is circularly polarized along the spin axis. The polarization becomes 
more elliptical as the angle from the spin axis is increased, becoming linearly polar- 
ized in the plane of the equator. Since the Faraday fading becomes the most severe 
for a linearly polarized wave, the observance of this severe fading indicates the station 
to satellite vector is normal to the spin axis. 

During the first two months, the orientation appeared to remain almost constant 
at the injection condition. A slight shift was indicated, but the limits of accuracy of 
the measurement were not sufficient to ascertain this. Thermal measurements taken 
throughout the first year of orbiting were in general agreement with the assumption of 
a space-fixed orientation, although again such measurements are not sensitive enough 
t-o detect a slight shift of the spin axis [Ref. 16]. 

It is certain that gravitational torques act on Explorer VII and it seems virtually 
impossible to construct a satellite without some magnetic moment unless particular 
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care is taken to prevent such a moment . Explorer VII does have a very large angular 
momentum resulting from its fast spin rate which greatly reduces the forced precession 
of the spin axis. 

The motion of the spin axis was integrated, assuming only gravitational torques, 
and resulted in the dashed curve in Figure 10. This is seen to result in a very slow net 
drift of the spin axis with time. The addition of a small magnetic moment with M 3 B^. = 
75 dyne-cm produces a motion indicated by the solid line. Again, it is seen that al- 
though an oscillation of several degrees is present, the net drift is quite small. This 
motion also appears to fit the observed orientation data, although the accuracy is not 
really sufficient to confirm it. 


SECTION VII. CONCLUSIONS 


It now seems clear that permanent magnetic and gravitational torques are the 
dominating effects which produce changes in a satellite orientation. Much the same 
conclusion was reached by Manger [Ref. 17]. for the Tiros satellites. In fact, the 
orientation of Tiros III was effectively controlled by applying currents in a loop to alter 
the magnetic moment. The orientation was predicted and the control currents were 
programmed prior to launch. According to Hetch and Manger [Ref. 18] , this was 
quite successful. Other satellites have been passively oriented along the magnetic 
field lines by employing permanent magnets and hysteresis damping rods [Ref. 19]. 

The magnetic effects in Explorer XI and Explorer IV were not as obvious as those 
in stable spinning satellites since it was not clear that a permanent magnetic moment 
projection on the axis of rotation existed. But now that the motions of these satellites 
have been successfully analyzed and understood, aheretofore disturbing problem has 
been resolved . 

It is now felt that sufficient knowledge of the rotational behavior of satellites exists 
to make it worthwhile to formulate the equations of motion about the center of mass 
with those for the orbit and solve the full six: degree of freedom equations of motion. 
Particular care should be taken to determine the moments of inertia and magnetic 
moments prior to flight. 

It was also shown that simplifying approximations may be employed that give a 
fairly good representation of the motion without integrating the orbital equations. The 
orbital information necessary to use these approximate equations is the semi-major 
axis, the inclination, and right ascension of the ascending node. These equations are 
simple enough to be integrated on a small computer such as a ReComp II or LGP-30 
or, if necessary, by hand. 
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FIGURE 10. COMPUTED ORIENTATION OF EXPLORER VII WITH AND WITHOUT 
AN ASSUMED MAGNETIC MOMENT 



Knowledge of the torque-producing effects is also of vital importance in active . 
attitude-controlled satellites since such information can be used to determine the op- 
timum saturation level for inertial devices and propellant storage requirements for 
mass expulsion devices. 
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